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Abstract 

We study the evolution of entangled coherent states of the two quantized elec- 
tromagnetic fields under dissipation. Characteristic time scales for the decay of 
the negativity are found in the case of large values of the phase space distance 
among the states of each mode. We also study how the entanglement emerges 
among the reservoirs. 
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1. Introduction 

Entanglement of two or more systems is an important concept playing a 
central role in quantum information science and quantum computation [1, 2]. It 
is an essential resource for teleportation [3] , quantum key distribution [4, 5] and 
controlled quantum logic [6], among others. The study of entangled states of 
elementary quantum systems has become a fruitful field of research [1, 2, 3, 4, 5, 
6, 7, 8, 9, 10, 11]. In particular, it has been paid a great attention on studying 
the behavior of quantum entanglement as a function of time when the system 
is affected by the environment [12, 13, 15, 16, 17, 18, 19, 20, 21, 22, 23]. 

The connection between decoherence and disentanglement has received con- 
siderable attention in the last years. Whereas decoherence of systems in con- 
tact with Markovian reservoirs decays asymptotically, entanglement may dis- 
appear suddenly. This has been recognized both in continuous variable sys- 
tems [12, 13, 14] and in finite dimensional systems [15, 16, 17, 18, 19, 20]. In 
particular, entanglement among coherent states is a subject that has received 
attention by a number of authors Ref. [24, 25, 26, 27, 28, 29, 30]. 
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The purpose of this paper is to study entanglement properties of superposi- 
tions of two-mode coherent states, in contact with a dissipative reservoir paying 
particular attention to the general features of the dynamics of entanglement. 
When the phase separation of the coherent states of each mode is large, we 
show, in a large variety of numerical examples, that the entanglement decays 
rapidly with a well defined characteristic time scale, which is much smaller than 
the time, r, when the entanglement completely disappears. This is in contrast 
to the case of two two-level systems, where the value of t does have a physical 
relevance. Analytical results are also found for certain classes of states showing 
the same dynamical behavior. By modeling the reservoirs through a set of linear 
oscillators we also discuss how entanglement emerges among the reservoirs as 
the system looses it. 

This work is organized as follows: In Sec. II we present the principal features 
of the entanglement dynamics under dissipation when considering as the initial 
state a general superposition of two coherent states for each mode. In Sec. Ill we 
study particular situations where we may easily describe analytically the time 
evolution of the negativity and the conditions for complete disentanglement. 
In Sec. IV we consider explicitly the reservoir degrees of freedom, in order to 
describe how the entanglement emerges into the two reservoirs. Finally, in Sec. 
V we present our concluding remarks. 

2. Disentanglement dynamics 

Consider the dynamics of two dissipative quantum modes, each affected by 
its own environment. Such a situation can be conveniently described at zero 
temperature by the master equation 

p^Y (^"'^"^ ~ ~ '°°^") + y (^^'^^^ " ^^^'^ " '°^^^) ' 

where a, and b, are the annihilation and creation operators of the two 
modes. 

Consider the class of initially entangled pure states of the form 

|*)= J2 ^Uc^^)Wk), (2) 

i,fc=l,2 

where ja^) and \a^) and l/Jj) ) ^-re coherent states associated with mode a 
(mode b) and {i,k = 1,2), are complex constants. An important feature of 
the solution of Eq. (1) for p{t) is that, for each t it can be written only in terms 
of the operators \ai{t) , /3m{t)) Pn{t)\- In fact, it can be shown that this 

solution is given by 

m= E C™-"(t)|a.(i)/3™W)(a,W/3nWI, (3) 

'i,j,m,?T,— 1,2 

where a,{t) = a^e-^^'/^, l3i{t) = PiC-^"'^/'^ , « = 1,2, and 
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CTfit) = cr^")*exp{-(l - e-r»*)(|a, - a,f + 2ilm(aia*))/2}x 
exp{-(l - e-i^"*)!!/?™ - Pn? + 2iluY{P^I3l))/2}. 



(4) 



The state above could be conveniently written in an orthogonal basis for each 
time t. For example, we could use the basis 



la(0> = »^.(lf.W>±|f|l|l^^«>). (5) 

with = 1/^/2(1 ± 1(6 and ^ = a,(i. 

Although our system is a continuous one, it evolves in such way that at any 
time t it can be described in the time dependent subspace spanned by the vectors 
\a±{t)) ® |/3±(<)). In this basis, we are allowed to treat the whole system as an 
effective two-qubit system [31] so that we can study entanglement properties 
through the concurrence [32] 

C{t) = max{0, (Ai - As - A3 - A4)}, (6) 

where Ai(t), i = 1, 2, 3, 4 are the square root of eigenvalues, in decreasing order, 
of the matrix M{t) — p{t)aiy®cT2yp{t)* cTiy®a2y, with all matrices been written 
in the time dependent 4-dimensional basis [Eq. (5)]. 

The coherence of a superposition of a pair of coherent states decays for short 
times with a time scale depending on the distance between the superposed states 
and almost disappears much before the characteristic time for energy dissipation. 
The problem we want to address here is how entanglement is lost for a pair of 
entangled coherent states and how it depends on the initial conditions. For any 
state of the form of Eq. (3) there is no simple analytical expression to describe 
how the entanglement evolves. In this section we take a numerical approach to 
get some insight in understanding the evolution of entanglement as a function 
of distances between the superposed states. 

For short times compared with l/rQ(l/rb) the matrix elements C"^", i ^ j, { 

C™'", m ^ n,) decay exponentially with a time constant Ta = [(Falai — a2|^)/2] 

( Tfc — [(Fftl/?! — /32|^)/2] ^ ). Ta (rb) are time scales associated to the loss of 
coherence among the states [ai) and |a2)(|/32) of mode a (b). Entangle- 

ment, as measured by the concurrence, Eq.6, is a sum of non linear functions 
of the matrix elements C™'"(i), and, depending on the initial conditions, may 
even vanishes at finite times although the matrix elements C™'"(t), only vanish 
at infinity. On physical grounds, we expect that entanglement evolution will 
depend mostly on the matrix elements C,™'"(i) with i ^ j,m n. Therefore for 
values of D"^ = {\ai — 02!^ + — /32|^)/2 ^ 1, and for times shorter than l/Ta 
and 1 /Ft, the concurrence should initially drop exponentially in a characteristic 
time r = l/[(ra|ai — 02!^ -l-Ffc|/3i — /32|^)/2]. We show examples of this behavior 
by plotting, in Fig. 1, \nC{t)/C{0) versus t/r for Ta = T^, for fixed values of 
= (\ai — Qfsl^ + |/3i — /32|^)/2 and for values of the constants chosen at 
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random. The curves show evidence that the entanglement initiahy does follow a 
single exponential behavior. After a while some of the curves drop more rapidly 
showing evidence of a non asymptotic decay of the concurrence, suggesting that 
finite time complete disentanglement is present for a variety of states. In the 
next section we study special cases for which one may obtain simple analytical 
results for the concurrence. 



3. Special cases 

In this section we consider certain classes of states where we are able to find 
simple analytical results for the concurrence, showing explicitly the behavior of 
disentanglement suggested in the last section. In all cases we take, for simplicity, 
Fa = Fh = r and the amplitudes ai(/3i) and a2(/32) of the states in Eq. (2) 
having the same phase. The density matrix in Eq. (3) can then be written at 
each time t in the orthogonal basis defined by 

\±,t)a = Na±im a,{t))± \ 

\±,t)b = m±it){\P,{t))±\Mt))), (7) 



where 



Na±it) = l/J2(l±e-3(l"i-"^l')'="") 



Nb±{t) = l/V2(l±e-5 



-|(l/3l-&P)e- 



(8) 



We also use, from now on the negativity, defined as smallest eigenvalue of the 
partial transposition matrix p^^ [33, 34], as a measure of entanglement. 
As a first example we consider the state 

I*) = + \ail32) + \a2fh) - |a2/32», (9) 

which is a very high entangled state whenever c?^ = |ai — a2p/2 ^ 1 and 
d| = \(3i — /32p/2 ^ 1. Within these conditions we obtain a simple form for the 
negativity: 

A_ = (1 - e(^"'"*-i"^" - e(^"'^'-i)^^ - e(^"'"*-i)'^=e(^"''''-i)4)/4, (10) 

so that entanglement decays with a superposition of exponentials within a char- 
acteristic time 1/ r(d^ + d|) as suggested in section II. Numerical results for 

the negativity are shown in Fig. 2, for F^ = Fh = F and = (d^ + d|)/2 
equals to 0.5,0.75, 1.5,2 and 5. As we observe in the figure, In A_(<)/A_(0) as 
a function of TD^t has a slope close to —1, even for as small as 2. 

Consider now that the system is initially prepared in entangled states of the 
form 

1*) =c+,+ |+,0)J + ,0), + c_,_|-,0)J-,0),, (11) 
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where, without loss of generahty, c+,+ and C-- are taken as real and positive 
constants. The density matrix written in the basis {| ±±, t) :=\ ±, t)a® \ ±, t)b} 
has the form of an X — matrix as previously considered for studying finite-time 
disentanglement [18, 19, 20]. In such case it can be shown that the concurrence 
is related to the negativity, A_, by the expression C{t) = max{0, — 2A_}. For 
simplicity, let us assume equal distances for each mode, |ai — = — /32|- 
The negativity can then be easily calculated and it is given by 



A_= l-e-° 



Al - A?e-^'(i-^"")l , (12) 



with Ai = c+.+A^a+(0)7Vfc+(0)+c_,_iVa-(0)A^fc-(0) and ^2 = c+,+iV„+(0)iVb+(0)- 

C^_A^a-(0)iVb_(0). 

Complete disentanglement occurs when the eigenvalue A_ becomes positive. 
From the expression above it is easy to see that this happens in a time t^^^ given 
by: 



The condition for having finite disentanglement time is that exp(— Z?^/2) < 
IA2I / |Ai| , which can be satisfied in two different cases. Firstly for c___// > 
C+.+ , with / = (1 — e^^ ^^)/(l + /^) for which there is always a finite 
disentanglement time t^^^ if c_^_ > c_|_.+ . Secondly when c_|-.+ > c^,^/ f, 

where there is t'^p finite only if c___ < c+_+/^. 

For large distances, Z? ^ 1, we see from Eq. (12) that the concurrence drops 
exponentially with a characteristic time (D^F)"-'^, which is much shorter than 
the dissipation time F""'^, in agreement with the results suggested in section II. 

Although the loss of entanglement may happen at a finite time t'^^\ the 
concurrence, for D > 1, is already small at a much earlier time. A similar result 
was found in Ref. [35] when studying GHZ N qubits states under different 
independent reservoirs. Consider, for example, the expression (12) for the case 
of the initial state given by Eq. (11). For ^ 1, Tt^^^ <C 1 we can easily show 
that the negativity drops significantly to low values much before the time 

when it completely disappears, independently of the values of c++,c and D. 

In Fig. 3 we show the behavior of the negativity for a fixed value of the ratio 

c++/c and for several values of D^, as a function of the normalized time 

r = t/tl^\ From that figure we see that the negativity drops to very low values 
in a time less than half the complete disentanglement time t^. 

For small distances, D 1, Eq. (12) can be approximated to 



A_~c2 e-^' 



1 - ^ - e-" 
c 



(14) 



The interesting about this regime is that clearly shows two different scales 

of disentanglement depending on the values of c+j^ and c In particular, 

for an initial maximally entangled state, that is, c-|_+ = c = l/\/2, A_ ~ 

— (1/2) exp (— 2Fi) the entanglement decays asymptotically with a rate 2r. For 
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c___ > c+__|_ entanglement vanish for time td = — F"-'^ ln(l — c+,+/c_,_). For 
c_^_ < c+^+ entanglement decays always asymptotically. So that this case be- 
haves in a similar way to the two-qubit case in [20] . 

Other situation of interest is the class of initial states having the form: 



1$) - c+,_|+, 0) 0), + c_,+ I-, 0) J+, 0) 



b ' 



(15) 



where c+^_, c__-|_ are taking as real constants. The dynamics of this state leads 
also to an X density matrix. In a similar way as for the previous case, we found 
that the time for disentanglement is given by: 



.(2) 



1 



L>2 



In 



(2) 

The only condition required for having finite is that: 
1 



(c+,-+C-,+) 



eM-^D') < |(c+, 



C-,+)/(c+,- +c_ 



(16) 



(17) 



4. Birth of entanglement among the reservoirs 

In the above sections, we have studied disentanglement dynamics for some 
special cases of entangled coherent states, paying attention only to the system 
of the modes a and b. However, when including degrees of freedom of both 
reservoirs, a deeper understanding of entanglement dynamics can be obtained. 
This can be carried out by describing the system-reservoir dynamics per each 
mode with a Hamiltonian of the form [37] 

N N 

H = huja^ a + h'^ojkclck +h'^gk{acl -t-c^a^). (18) 

k=l fc=l 

In such case, it is not difficult to show that when the mode a is initially in a coher- 
ent state and the reservoir is in the collective vacuum state [0)^ = IlfcLi |Ofc)ri 
the mode a and its reservoir evolve into a product of coherent states of the 
form ] at)\at)r —\ oiP{t))Y\^^i\fk{t)ct)r- In the Born-Markov approximation 
F{t) = e-r»*/2 and Efe fk{t)^ = 1 " e"^"*. 

To illustrate how entanglement emerges into the reservoirs, let us consider 
the initial state of Eq. (11), but now including the states for both reservoirs 

]*) = (c+.+ 1+, 0)^ 1+, 0), + c_,_ 1-, 0), 1-, 0), ) 10),„ 10),, , (19) 

where c+^+ and c_ _ are real and positive constants. Note that a{b) and 7'a(f)) 
stand for mode a{b) and reservoir associated to mode a(b), respectively. Since 
the evolution of each mode and its reservoir may be treated independently, the 
complete density matrix at time t can be obtained through the evolution of the 
states l±,0)alO),„ = A^a±(0)(la)±l-a))lO),„ and |±, 0)6l0),, = iVb± (0) (]/?) ± 
|-/3»|0)... 
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Let c and Tc stand for either mode a or 6 and for their own reservoirs. It is 
possible to show that the initial states |+,0)c|0)r^ and | — ,0)c|0)r^ evolve at the 
time t to the states 

(20) 

respectively. The states |-E±.t)rc in the equations above are defined as |i?±,t)r„ = 
(|at)r„ ± I - dt)r^)/Ba±it) and \E±,t)n = {Wt)r^ ± | - /3t),.J/B6±(t) and fo°rm 
a two-dimensional orthogonal basis for each reservoir, and normalization terms 
read 



Ba±{t) = ./2(l±e-2|"P(l-e-ra*)) 



Bb±it) ^ ^2{l±e-mHi-e-^^')), (21) 

Let us now study the joint evolution of the two modes coupled to their 
reservoirs. In the above description, the initial state (19) evolves to 

l*)t = C++\4>+^t)ara\4'+,t)bri, + C ara\4'- ,t) brt ■ (22) 

This is a pure state representing the whole system. As t increases, there is a 
transfer of the correlationsof the system, both classical and quantum, to the 
reservoirs. Notice that there is no memory effects involved in this transfer. 
Information about the reservoir-reservoir dynamics are obtained by tracing out 
cavities degrees of freedom, as well as the cavity modes evolution are obtained 
by tracing out thestates of the reservoirs. 

As before, let us consider equal distances for each mode, = \ai — = 
\Pi — /32p- The dynamics given by Eq. (22) leads us to an X — matrix for the 
reservoir-reservoir subsystem, which is complementary to the density matrix of 
the subsystem of the two modes a and b [37] . As a consequence, we can study 
entanglement dynamics associated to external degrees of freedom. 

When studying entanglement through the partial transposition matrix, we 
obtain a complementary result of (12) giving us information about the birth of 
entanglement in the reservoir-reservoir subsystem. In this case, we obtain the 
following negative eigenvalue 



Ai - Aie-^ ^ , (23) 



from which we can derive the time for birth of entanglement in reservoirs, that 
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Note that Eq. (23) can be obtained from Eq. (12) by exchanging the time- 
dependent functions, e~'"* and 1 — e~^*. Also 

e-^*" + e"^*" = 1 (25) 

This result immediately shows that if the entanglement among the two modes 
persists all the time, the entanglement among the reservoirs starts growing at 
t = 0. Eqs. (13) and (24) shows that in this case the entanglement among the 
reservoirs never dies and its negativity reaches at infinity the initial value of 
the negativity among the two modes a and b. Eq. (25), together with Eqs. (13) 
and (24) also shows that the time when entanglement starts among the two 
reservoirs can be either smaller, greater or equal to the time when entanglement 
among the two modes completely disappears. This result is analogous to the 
one obtained in Ref. [37]. 

Figure 4 shows the entanglement evolution when tf, > td, considering two 
distances in phase space = 0.4 (Fig. 4(a)) and D'^ = 16 (Fig. 4(b)). Here 
there is a time window where no entanglement exist either among the two modes 
or among the two reservoirs. However, we observe that entanglement does exist 
in others bipartite partitions of the whole system, since each mode entangles 
with the reservoir of the other mode. When plotting the concurrences Cab, Cr^rj,, 
Cara S'lid Car^, ^c noticc that in the region where there is no entanglement 
either among the modes a and b or among the two reservoirs, that is. Cab = 
Cr^rb = 0, entanglement between a mode and its corresponding reservoir reaches 
its maximum value. Although this behavior is presented in both cases, the 
correlations are sensitive to any change of parameter , as shown in Fig. 4(b) 
where the concurrence between mode a and reservoir rb is very small. Moreover 
for small distances (D^ = |2Q:p), which also correspond to a small average 
number of photons, we reproduce the entanglement dynamics associated to the 
two-qubit case studied in Ref. [37]. Note that since we are considering equal 
distances and decay rates, both mode-reservoir subsystems evolves identically 
implying that the concurrences Car^ = Cbn, and Car^ = Cbr^ ■ 

A particular result happens when IA1/A2I — that is c-|_,+/c___ = 

(1 + e^^ /^)/(l — e^^ /^Y. In this case, the times when the modes completely 
disentangles and the time when the reservoirs starts to be entangled are equal 
and occurs at h = td = ln(2)/r, a result which in independent of the distances 
in phase space for both modes. 

5. Conclusion 

In conclusion, we have studied the evolution, under dissipation, of the entan- 
glement among two modes of the electromagnetic field for certain class of ini- 
tially entangled coherent states. At zero temperature, the density matrix can al- 
ways be written, for each t, in a finite orthogonal basis, which allow us to describe 
the entanglement evolution as that happening in a discrete system. Asymptotic 
entanglement decay as well as finite time disentanglement arises depending on 
the initial conditions for the bipartite system and on the distances in phase 
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space among the components of each mode. Typicahy, the robustness of entan- 
glement decays in a time scale proportional to 2/iTa\oii — Q;2p + rb|/3i — /32p), 
the inverse of a weighted average of the distances in phase space for each mode, 
whenever these distances are much greater than one. Both numerical and ana- 
lytical results are presented. 

We have also presented an investigation concerning the entanglement gener- 
ation among the reservoirs. By modeling the reservoir and its interaction with 
the system, it is possible to study death of entanglement among both modes of 
the field and the birth of entanglement in reservoir-reservoir subsystem, respec- 
tively. In addition to the dependence on initial amplitudes, the characteristic 
times for death and birth of entanglement depends on the distances among co- 
herent states. An interesting result shows that for specific conditions in the 
initial amplitiides for a given distance among coherent states, both birth and 
death of entanglement at finite times occur simultaneously at a time which does 
not depend on the distances of each mode. 
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Figure 1: lnC(i)/C(0) as a function of t/r, for = 0.75 (line and circles), = 2 ( dashed 
line) and D-^ = 5 (solid line ). For each value of , three sets of coefficients c'y that define 
the initial state (2) were chosen at random. Fa = F;, = F 




Figure 2: lnA(t)/A(0) as a function of the renormalized time 

t [(Fa|«i - a2p + F^l/?! - /32p)/2] , for the state given by Eq. 9. = 0.5 (dashed 

line), = 0.75 (dotted line), = 1.5 (dot-dashed line), = 2 (circles), = 5 (solid 
line). Fa = Fi, = F 
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Figure 3: Ratio of the negativity at time t and at the initial time as a function of the normalized 

time r = t/t^^\ for fixed |c /c++| = y/2 and distances = 0.5 (dashed line), = 2 

(dotted line) and D-^ = 7.5 (dot-dashed line). It is also shown, for comparison, the line 1 — r. 
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Figure 4: Evolution of concurrence for different bipartite partitions. Cc^c^ (solid line), Cr^ri, 
(dashed line), Cc„r„ (dotted line), and Cc^r^ (dot-dashed line), (a) considering = 0.4 and 
(b) for D-^ = 16. In both simulations we considered equal decay rates Fa = F;, = F, and 



amplitudes are c+,+ = 0.4, c___ = 
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